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Abstract 

In this paper we apply the differential inequality technique of Payne et. al HU] 
to show that a reaction-diffusion system admits blow-up solutions, and to determine 
an upper bound for the blow-up time. For a particular nonlinearity, a lower bound 
on the blow-up time, when blow-up does occur, is also given. 


1 Introduction 


Since the pioneering work of Fujita ([I]) on the blow-up of solutions of nonlinear diffusion 
equations, there has been considerable interest in the study of such solutions for nonlinear 
parabolic equations. Blow-up phenomena can be observed in nature and is important in 
various disciplines such as biology, chemistry and physics. More detailed results related 
to the blow-up of solutions can be seen, for example, in the monograph [T^l [131 US], the 
surveys [Diaiisiiii] and bibliographies cited therein. 

Recently, Payne and Schaefer (illZllEllgl HDl E]) have applied the energy method to 
derive a differential inequality for certain integrals corresponding to the blow-up solutions 
of parabolic equations. By means of this method, they find lower bounds on the blow-up 
time in certain nonlinear parabolic problems. In the present paper, we find that the above- 
mentioned method can be applied to reaction-diffusion systems to obtain generalizations 
of results in [10]. More precisely, it can be shown that Theorem 2.1 and Theorem 3.1 in 
[To] can be extended to Theorem 11.11 and Theorem 13.11 respectively in this paper. 

Compared with cases of single equations, little work appears to have been devoted 
to systems. Motivated by Theorem 2.1 in HD], we show that the result in Theorem 2.1 
can be extended to the gradient system case. To be specific, we consider the following 
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initial-boundary value problem for the reaction-diffusion system under a Robin boundary 
condition, i.e. 

Ut = Au + fi{u,v), in fix( 0 ,cxo), 

Vt = Av + f 2 {u,v), in f 2 x( 0 ,cxo), 

^ ( 1 - 1 ) 

gM+yiti = 0, ^ + -f 2 V = 0, on 51dx(0,CX)), 

^w(a;, 0 ) = 5 fi(x), v{x,0) = g 2 {x) in Q, 

where {u, v) = {u{x, t),v{x, t))] fl E > 2 is a bounded domain with smooth dfl] 71 

and 72 are positive constants; u is the unit outward normal on dfl; and gi{x) and g 2 {x) 
are nonnegative functions which do not completely vanish. Suppose that the nonlinearity 
fi{u,v) and f 2 {u,v) in fll.ll) satisfy 

(HI) ufi{u,v) + vf 2 {u,v) > 2(1 -F a) F{u,v), 

where F{u,v) is a solution to duF{u,v) = fi{u,v), dyF{u,v) = f 2 {u,v), and a > 0 is a 
constant. In addition to the non-negativity of the initial conditions gi{x) and g 2 {x), we 
impose the additional conditions: 

(H2) 2/^ F(^i(x),^2(a:))da; > 71 Jd^gfds + J^ \Vgi\^dx; 

(H3) 2 F{gi{x),g 2 {x)) dx > 72 fg^ gl ds + | V^ 2 p dx. 

Under these conditions we have: 


Theorem 1.1. Suppose that {u{x,t),v{x,t)) is a pair of solution to fll.ip . //(hTl) (H3) 

are satisfied, then at least one ofu{x,t) and v{x,t) blows up in finite time t*, where t* is 
bounded above by 02 . 101 ) . 


The physical meaning of the Robin boundary conditions can be explained as follows. 
Suppose that u and v represent temperature, and are governed by the equations in the 
problem 01.ip . Then the Robin boundary conditions mean that the heat flux ^ and ^ 
on the boundary of 12 are proportional to the temperature u and v on the boundary of 
12, respectively. Since 71 and 72 are positive constants, it follows that the larger the heat 
flux is, the smaller the temperature is. We note that, from the biological point of view, 
the temperature and the heat flux can be substituted respectively to population density 
and population flux. 

In other words, the larger the population flux is, the smaller the population density 
is. As a consequence, when the population flux on (212 is large, the population density 
on (212 is small. The low density of u and v on (212 then may result in the blow-up of u 
or V since the large flux flows into 12 but on the boundary of 12 , the density of u and v 
are restricted to be small. Therefore, u and v are may be forced to aggregate together so 
that blow-up occurs. 

The remainder of this paper is organized as follows. In Section |2l we give the proof 
of Theorem o This theorem asserts that for certain initial conditions, the solutions of 
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fll.lj) blow up in finite time, with an upper bound given by fl2.10j) . Section [3] is devoted to 
determining a lower bound on the blow-up time when blow-up does occur. In addition, 
Theorem 13.21 provides a cooperative system considered in [2] a lower bounded for blow-up 
time when blow-up does occur. Finally, we conclude the present paper with some remarks 
in Section m 


2 Proof of Theorem 11.11 


In this section. Theorem 11.11 is proven. 


Proof. First of all, we define 




( 2 , 1 ) 


By means of integration by parts and {HI), we arrive at 

E'{t)=2 I u{Au +fi{u,v)) + V {Av + f 2 {u,v)) dx 

Jn 

= —271 / vfds — 2 / \Vu\^dx + 2 / ufi{u,v)dx 
Jan Jn Jn 

— 272 / v^ds — 2 / \Vv\^dx + 2 / vf 2 {u,v)dx 
Jan Jn Jn 

> — 2 (1-|-a) (71 / ds + / I Vm|^ dx) — 2 (1 -|-ck) (72 / v‘^ds+ |Vn|^dx) 
Jan Jn Jan Jn 


-I-4 (1-1-a) / F{u,v)dx. 

Jn 


( 2 . 2 ) 


Letting 


J(t) = — 2 (1-|-a) (71 / u^dsP I Vnp dx) — 2 (1-|-a) (72 / v‘^ds+ |Vnpdx) 


Jan 

-|-4 (1-|-ck) / F{u,v)dx, 
Jn 


'an 


(2.3) 
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we calculate the derivative of J{t) to obtain 


J'(t) = — 4 (1 + a)( 7 i / uutds+ / Vu-Vutdx— / fi{u,v)utdx) 


I on 


4(1 +a )(72 / vvtds+ / Vv-Vvtdx 
Jan Jn 

du 


f2{u,v)vtdx) 


= - 4 (l + a)( 7 i / uutds+ I -Utds- Ut{Au + fi{u,v)) 


Jan Jan 

- 4(1 +a )(72 / vvtds + 

Jan Jan 

>4 (1 + a) / {ul + v 1 ) dx, 


dv 

du 


vtds- / Vt{Av + f 2 {u,v)) 


( 2 . 4 ) 


by virtue of integration by parts and the Robin boundary conditions in fll.ip . Due to 
{Jd2) and {H3), we have J(0) > 0. Then J(0) >0 and J'{t) > 0, for f > 0, yield J{t) >0 
for t > 0. Because 


E'{t) = 2 / {uut + vvt)dx 

Jn 

we have by applying Cauchy-Schwartz inequality 

=4 [ uutdx'^ + (^ I vvtdx'^ +; 


(2.5) 


u Ut dx 


V Vt dx 


<4 


dx uidx+ dx / v: dx 


w dx / ujdx dx / v: dx 


m 

<4 / (n^ + dx I {lif + v‘f) dx 


n 

2 I ..,25 


<- 


1 + CK 




(2.6) 


The last inequality holds since \/^ < for a, 6 > 0. From the dehnition of J(t), 
E'{t) > J(t), and consequently we can replace E'{t) by J{t) in fl2.6p to obtain 


J{t) E'{t) < 


1 -f* Q 


m j\t) 


or 




E(t) - J(() ^ 

Following the same arguments in [10], we obtain that +(t) satishes the inequality 

1 1 


< 


(E(i))“ - (£(0))‘ 


— a M t, 


(2.7) 

(2.8) 

(2.9) 
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where M = (^(o))i+a ■ Since the last inequality cannot be true for all t > 0, we infer that 
at least one of u{x,t) and v{x,t) blows up in hnite time, t*, where t* is bounded above 
by 


r < 


1 

aM{E{0))^' 


( 2 . 10 ) 

□ 


We note that Theorem 11.11 remains true if the Robin boundary conditions are replaced 
by Neumann boundary conditions (i.e. 7i = 72 = 0). To illustrate the results in The¬ 
orem [m we give an example of the nonlinearity in fll.ip . Take F{u,v) = then 

fi = duF{u,v) = 2uv^, f2 = dyF{u,v) = 

u fi{u,v) + V f 2 {u,v) = v^, (2.11) 

and 

2{l + a)F{u,v) = 2{l + a)u^v^. (2.12) 

Clearly, (HI) is fulhlled when 0 < a < |. 

Remark 2.1. Equivalently, {HI) can be rewritten as 

u duF{u, v) + V dyF{u, v) > 2(1 -|- a) F{u, v). (2-13) 

For equality, that is 

u duF{u, v) + V dyF{u, v) = 2(1 -|- a) F{u, v), (2-14) 

which can be solved by the method of characteristics to give 

F(ri, n) = /i^—^ , (2.15) 

where h = h{w) is an arbitrary smooth function and c is an arbitrary constant. This 
solution is useful in looking for nonlinearities /i and /2 in fll.ip which satisfy {HI). 


3 Lower bound for blow-up time 


In this section, a lower bound on the blow-up time is obtained when blow-up does occur. 
In particular, (2.16) in [S] plays an essential role in proving the following 


Theorem 3.1. Let {u{x,t),v{x,t)) be a pair of nonnegative solutions to fll.ip and at 
least one ofu{x,t) and v{x,t) blow up in finite time, t = t*. Suppose that (^41) (kl3) 

below are also satisfied: 
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(Al) Q E is a hounded smooth convex domain; 
(A2) fi{u, v) < ki for ki,u,v > 0 and p > 1; 

(A3) /2(w, v) < k 2 for k 2 , u,v > 0 and p> 1. 

Then t* is bounded below by fl3.12p . 

Proof. First let us define the auxiliary function 


^(^) = / + v'^^) dx. 


(3.1) 


On applying Green’s first identity and the equality = p , we arrive 

at 

£'{t) =2p f (Au + fi{u,v)) + {Av + f 2 {u,v)) dx 

Jn 

= — 2 p 7 i f ds — 2p {2p — 1) f \Vu\‘^ dx + 2p f fi{u,v) dx 

Jan Jn Jn 

— 2p72 f ds — 2p {2p — 1) f I Vnp dx + 2 p f f 2 {u,v) dx 

Jan Jn Jn 

< — 2 {f2p — 1) p~^ I |Vn^pdx + 2p j fi{u,v) dx 

Jn J n 

— 2{2p—l)p~^ f \Vv^\‘^dx + 2p f f 2 {u,v) dx 

Jn Jn 

< — 2{2p—l)p~^ / dx + 2p A:i / dx 


— 2{2p—l)p ^ / iVn^pdx + 2pA:2 / n^^dx 

J n J n 


(3.2) 


since (242) and (^43) hold. Now our strategy is to relate dx in terms of S{t) and 

iVn^pdx. To this end, the integral inequality (see (2.16) in [S]) is used: 


[ dx < —3- < 7^ [ dx + (- + 1 

'o 34 1 2p \p 


u 


dx 


iVn^pdx ) , (3.3) 


where for some origin inside G, 


p = min Xj z/j, df = max Xj Xj, 
an Q 


(3.4) 
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for Vi the i—th component of the unit outward normal to dVL. Thus, 


£'{t) < — 2 { 2 p—l)p / \VuFy dx— 2 { 2 p—l)p / 


34 1 2 p Jq \p 



d 


(“ + 1 

34 \ 2 pJn \P 




\VuP\‘^dx 
\VvP\^dx 


(3.5) 


By virtue of the elementary inequalities (0 + 6)2 <22 (02 + 62) and 04 64 <^a + |6, we 
obtain 

£^'(f) < — 2 (2p — l)p“^ [ |Vu^|^ dx — 2 (2p — 1) [ \Vv^\‘^ dx 

dx 



dxj + ^- + 1 


dxY + (- + l"" 






P 



v^P dx 



\VuP\^dx 


\VvP\^dx 


< — 2 { 2 p—l)p ^ / |Vu^f dx — 2 (2p — 1)p ^ / iVyPfdx 

Ju Jq 



u^P dx 

/3r^ 




iVu^pdx 



v'^P dxj 


3 


v^Pdx'] 





(3.6) 


where and (32 are positive constants satisfying 


, , 1 S^pki f d ^ , ,1 3 ipk 2 f d \ ^ ^ 

-2 2p-l p-i + —^ - + 1 / 3 i< 0 , -2 2p-l p-^ + —^ - + 1 /?i<0. 

22 \p ) 22 VP / 

(3.7) 


7 
















Accordingly, 


, 2pki 1 

S'{t) <-^22 


3t 



34 






(3.8) 


3 3 

<34 pk p~^ 



+ 


pk [ d 


2332 


P 


+ 1 /3 


-3 


dx 


In 


(3.9) 


where k = max(/ci, /c 2 ) and f5 = niin(/5i,/ 32 ). As a resnlt, we obtain 


S'{t) < K,SHt)+K2S^{t), 


where 


Integrating yields 


Ki = 3ipkp 2 ^ K2 = ^^(- + A /3 

2234 \p J 


t* > 


dC 


Isio) Ki^2+K2e 

This completes the proof of the theorem. 


(3.10) 

(3.11) 

(3.12) 

□ 


We remark that Theorem 13.11 remains trne if the Robin bonndary conditions are 
replaced by Nenmann bonndary conditions. From the proof of Theorem 13.11 it is readily 
seen that the following resnlt is trne nnder an assnmption which is weaker than (A 2 ) and 
(A3). 


Theorem 3.2. Let {u{x,t),v{x,t)) be a pair of nonnegative solution of fll.ip and that at 
least one ofu{x,t) and v{x,t) blow up in finite time, t = t*. Suppose that (Al)' (A 2 )' 

below are satisfied: 

(Al)’ is a bounded smooth convex domain; 

(A2) ’ For ki, /c 2 , m, n > 0 and p>l, 

/ i ( m , v) + f 2 {u, v) < ki + k 2 (3.13) 

Then t* is bounded below by fl3.12p . 


Now we are in the position to apply the following results to give an example of the 
above theorem. 


















Theorem 3.3 (From [2]). Consider the following reaction-diffusion system with absorp¬ 
tion: 

Ut = Au-\-— au^, in f2x(0,cxo), 

Vt = Av + — bv^, in fix(0,cxo), 

< (3-14) 

u = 0, V = 0, on dflx{0,oo), 

^u{x,0) = gi{x), v{x,0) = g 2 {x) in n, 

where G M"'; a, b, p, q, r and s are positive constants; gi and g 2 are nonnegative 
functions. 

(i) Ifpq > max(r, 1) max(s, 1), then there exist solutions of fl3.14p which blow up in 
finite time. 

(a) Ifpq < max(r, 1) max(s, 1), then all solutions of fl3.14p are global. Moreover, if 
f'C > 1 (hence pq < r s), they are uniformly bounded. 

(Hi) If pq = max(r, 1) max(s, 1), then 

(a) ifr,s>l and a and b are sufficiently small, then there exist solutions of 03.14p 
which blow up in finite time; 

(b) ifr,s> 1, and a'^ 6'’ > 1 (equivalently, a^lf >1), then all solutions are global 
and uniformly bounded; 

(c) if r ors<l, then all solutions are global (possibly bounded). 

For the nonlinearity fi = — a and fi = — bv^ (by choosing p = q = r = s = 3 

in Theorem 13.31) . we let p = 2 in Theorem 13.21 so that 

fi{u,v) + f 2 {u,v) = - au^ - bv^, (3.15) 

while 

kiU^P + k2V^P = + k2V^. (3.16) 

By choosing ki, k 2 > 1, it is easy to see that 03.131) antomatically holds for any a,b > 0. 
We apply Theorem 13.31 (Hi) — (a) to conclude that 03.141) with p = q = r = s = 3 and a, b 
sufficiently small admits a blow-up solution. Now Theorem 13.21 provides a lower bounded 
for blow-up time t = t*. 


4 Concluding Remarks 


In Theorem 11.11 and Theorem 13.11 the estimates of the blow-up time for the blow-up 
solutions of 01.ip are given by means of the differential inequality technique adapted 
from mil- In the present paper, we have shown that this technique can also be applied 
to certain reaction-diffusion systems. From the proofs of Theorem 11.11 and Theorem 13.11 
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however, we cannot determine whether one of u and v blows up or if both u and v blow 
up. This problem is left for future work. 
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